WASHINGTON UNIVERSITY 


MEMORANDUM NO, 72/4 


OCTOBER., 1972 


COMPUTER PROGRAMS FOR PLOTTING SPOT-BEAM 
COVERAGES FROM AN EARTH-SYNCHRONOUS SATELLITE 
AND EARTH-STATION ANTENNA ELEVATION ANGLE CONTOURS 


PIECES SUBJECT TO ClMcIi 




# 


THOMAS W, STAGL 
JAI P, SINGH 


*JS 

' v*' 


% 

\-U\ 

U y 




Reproduced by 

NATIONAL TECHNICAL 
INFORMATION SERVICE 

US Department of Commerce 
Springfield, VA. 22151 



WASHINGTON UNIVERSITY 


ST. LOUIS / MISSOURI 


6 3 13 0 



PROGRAM ON APPLICATION OF COMMUNICATIONS SATELLITES 
TO EDUCATIONAL DEVELOPMENT 


Center for Development Technology 
Washington University 


Memorandum No. 72/4 


October, 1972 


COMPUTER PROGRAMS FOR PLOTTING SPOT-BEAM 
COVERAGES FROM AN EARTH-SYNCHRONOUS SATELLITE 
AND EARTH-STATION ANTENNA ELEVATION ANGLE CONTOURS 


THOMAS W. STAGL 
JAI P, SINGH 


This research is supported by the National Aeronautics and Space 
Administration under Grant No. Y/NGl-26-008-054. This memorandum 
is primarily for internal distribution and does not necessarily 
represent the views of either the research team as a whole 
or the National Aeronautics and Space Administration. 


V 



TABLE OF CONTENTS 


NO. Pa 9 e 

I. Introduction. 1 

II. Data Tape 2 

III. Plotting Programs 3 

IV. Antenna Coverage Patterns Using Data Tape . 4 

V. Polar Perspective 5 

VI. Footprints of ATS-F S-Band ETV Off-Axis 

Dual Feeds 7 

VII. References. 8 

Appendix A: WORLDMAP Program 

Appendix B: MINMAP Program 

Appendix C: ANTOVLY Program 

Appendix D: PERSPECT Program 

Appendix E: ATS-F Program 

Appendix F: Input Formats 



LIST OF FIGURES 


No. Page 

1. Flowchart of WORLDMAP Program 3a 

2. Flowchart of MINMAP Program 4a 

3. A Typical WORLDMAP Plot with an Antenna 

Coverage Overlay 4b 

4. A Typical MINMAP Plot with an Antenna 

Coverage Overlay 5a 

5. Polar Perspective Geometry 5b 


6. A Polar Perspective of the Earth From a 

Satellite at 120°W Longitude and Earth- 
Station Antenna Elevation Contours 6a 

7. Flowchart for the Polar PERSPECTive Program. 6b 

8. ATS-F S-Band ETV Footprints 7a 





I. INTRODUCTION 

This memorandum* contains a description and listings of 
computer programs for plotting geographical and political 
features of the world or a specified portion of it, for plot- 
ting spot-beam coverages from an earth-synchronous satellite 
over the computer generated mass, and for plotting polar 
perspective views of the earth and earth-station antenna ele- 
vation contours for a given satellite location. The programs 
have been prepared in connection with a project on Application 
of Communication Satellites to Educational Development. This 
report is intended to be a continuation of the work earlier 
reported by Stagl and Singh. 

A data tape was obtained from the Rand Corporation con- 
taining a digitized map of the world. Several programs have 
been written utilizing this data tape. The first program, 
called WORLDMAP, simply plots the tape as a map with axes and 
grid lines. The second program, called MINMAP, plots a speci- 
fied part of the total map so that larger scale plots of a 
small area of interest can be obtained. The program described 
in Memorandum 72/3 (1) is used as an overlay on these two plots 
This program plots the locus of intersection of quadric cones 
(narrow-beams from satellites) and a sphere (the earth) . 

A third program utilizing this data tape, called PERSPECT 
plots a polar perspective view of the earth and earth-station 
antenna elevation angle contours. This is a view of the earth 
as seen from a satellite in the geosynchronous orbit. Using 
a plot of this kind, one can make an overlay of the cross- 
sectional shape of the desired antenna beam and this overlay 
will be valid over the entire plot of the earth. This program 
is intended to facilitate determination of the area coverage 
by satellite-borne shaped beam antennas. 


*The authors wish to acknowledge the assistance of Mr. Neil 
Ostrander of the Rand Corporation, Santa Monica in acquiring 
the computer tape containing the digitized map of the world. 
They also wish to thank Mrs. Emily S. Pearce and Ms. Donna 
Barnes for typing the various drafts and the final version 
of the manuscript. 
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Also included in this memorandum is a modified antenna 
coverage program to be used to plot the footprints produced 
by the two off-axis S-band ETV feeds on Applications Technology 
Satellite-F (ATS-F) . This program can also be used as an 
overlay for WORLDMAP and MINMAP . 


II. DATA TAPE 

The point coordinates are stored in strings; latitude, 
longitude, latitude, longitude. . . . 0.0, 0.0. The range 
in latitude is from -it/2 to ir/2 (radian measure) , and the 
longitude range is from -it to it. The tape is formatted into 
24 blocks. The first block contains 23 integers which specify 
the number of coordinate pairs in each of the following 23 
blocks. Each of the following blocks contain a number of co- 
ordinate strings. The strings are of variable length and all 
strings end with a (0.0, 0.0) coordinate pair. Political 
boundaries are distinguishable from geographical boundaries 
in the following way. Signal coordinate pairs appear at the 
beginning of groups of strings. If the first coordinate pair 
in a string is (4.0, 0.0) then that string and all following 
strings are geographical boundaries until a string whose first 
coordinate pair is (8.0, 0.0) is encountered. This indicates 
that that string and all following strings are political 
boundaries until a string whose first coordinate pair is 
(4. 0,0.0) is encountered. 

The original tape obtained from the Rand Corporation was 
an 8 track 1600 bpi tape. Since the IBM 360/50 installation 
at Washington University Computing Facilities has no facili- 
ties for dealing with 1600 bpi tape, the data tape had to be 
copied onto an 800 bpi (bits per inch) tape which is kept in 
the Washington University Computing Facilities. A card backup 
was also obtained in case the tape is lost or inadvertently 
modified. 
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III. PLOTTING PROGRAMS 

The first program written to plot the data tape, called 
WORLDMAP, does not take advantage of the ability to distinguish 
between geographical and political boundaries. This is because 
the added programming complexity needed to distinguish political 
boundaries from geographical boundaries on the Calcomp plotter 
is not warranted by the use of the plot. 

Basically the program reads the coordinate strings from 
the tape and connects the points with straight line segments. 

The only problem is that the boundary of Siberia "wraps around" 
the end of the map and those strings that "wrap around" must 
be broken in order to avoid "retrace" lines across the plot. 
This is done by checking consecutive points in a string for a 
separation of over 2 radians. If the separation is greater 
than 2 radians the points are assumed to "wrap around" the ends 
of the map. In this case the string is broken there and plotted 
The remainder of that string is treated as a new string. 

The input to the program, FACT, is a magnifying factor 
for the map. With an input of FACT =1.0 a map 50 inches by 
26 inches will be drawn. Any number greater or less than 1.0 
can be used but one should consider the size of the output 
plot and the size of the paper available. Figure 1 is a 
flowchart of the WORLDMAP program. A listing of the program 
is included in Appendix A. 

The second program written for this data tape, called 
MINMAP, plots only a portion of the total map. The inputs 
to the program are the upper and lower limits of the longitude 
and latitude and the magnifying factor as for the WORLDMAP 
program. For this program, all points of each string are 
checked against these limits and either accepted or rejected 
depending on whether or not they are within these boundaries. 

The size of this plot is somewhat dependent on the input 
boundaries. The vertical axis is set at 10 inches. Using 
this fact and the upper and lower limits on the latitude, a 
scaling factor is calculated. This scaling factor is then 
















used, along with the upper and lower limits on the longitude, 
to calculate the horizontal axis length. As in the WORLDMAP 
program, the magnifying factor input can be used to change 
the size of the plot. 

The computer prints out the values of the first longitude 
value on the horizontal axis (FVX) , the first value on the 
latitude axis (FVY) , and the scaling increment (DA) in degrees 
per inch. These values will be needed for proper scaling of 
the antenna coverage program. 

A listing of the program is included in Appendix B. 

Figure 2 is a flowchart of the MINMAP program. 


IV, ANTENNA COVERAGE PATTERNS USING DATA TAPE 

The antenna coverage patterns program of Memorandum 72/3 
was modified for use with the WORLDMAP program. The modifi- 
cation, called ANTOVLY, includes an elimination of the portion 
of the program that scales the coordinates and draws the 
axes, as well as a redefinition of the origin and scaling 
factor to those of the WORLDMAP program. The modified pro- 
gram is listed in Appendix C. When the programs are run, 
the map is first plotted. Then the antenna coverage patterns 
are plotted over the map. The resultant plot shows explicitly 
the areas covered by the antenna beams. 

The inputs to the program are the first longitude value 
on the horizontal axis, FX, the first latitude value on the 
vertical axis, FY, the horizontal scaling increment, DX, and 
the vertical scaling increment, DY. The magnifying factor, 
FACT, is also input. 

When used as an overlay to WORLDMAP, the inputs should 
be: FX = -180., FY = -90., DX = 7.2, DY = 6.923, for matching 

the size of ANTOVLY plot with the WORLDMAP. The magnifying 
factor should be the same as that used when WORLDMAP was run. 

A sample of this type of output is shown in Figure 3. 




Figure 2. Flowchart of MINMAP Program 
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The ANTOVLY program can also be used with the MINMAP 
program. When used in this way, the inputs should be: 

FX = FVX (output from MINMAP) , FY = FVY (output from MINMAP) 
and DX = DY = DA (output from MINMAP) . The magnifying 
factor should be the same as that used when MINMAP was run. 
A sample of this output is shown in Figure 4. 


V, POLAR PERSPECTIVE 

The third program to use the data tape, called PERSPECT, 
was written to provide a view of the earth as seen from a 
satellite in an earth-synchronous orbit and to plot earth- 
station antenna elevation angle contours. Using a plot of 
this kind, one can make an overlay of the cross-sectional 
shape of the desired antenna beam and this overlay would be 
valid over the entire plot. The angular beam width would be 
proportional to the linear scale of the plot. 

To realize this type of plot, a projection is first made 
onto a sphere centered at the satellite as shown in Figure 5. 
The angle of the horizon, RMAX, is first calculated. A co- 
ordinate string is then read from the tape and the angle R 
is calculated for the first point according to the following 
formula : 

R = Cos” 1 [Sin (LTSS) Sin(YL) + Cos (LTSS) Cos (YL) Cos(XL)] 

( 1 ) 

where LTSS is the subsatellite latitude, XL and YL are the 
longitude relative to the subsatellite longitude and the 
latitude, respectively, of the point in questions. If 
| RMAX-R | < 30°, that is, if the first point of that string 
is within 30°, of the horizon, then the angle R is calculated 
for each point of the string according to equation (1) and it 
is checked for visibility. If (RMAX-Rj >_ 30°, that is, if 
the first point of that string is greater than 30° from the 
horizon, it is checked for visibility. If that point is 



LONGITUDE 



Figure 4. A Typical MINMAP Plot With An Antenna Coverage Overlay 





Spherical 



Figure 5. POLAR PERSPECTIVE GEOMETRY 
(not to scale) 
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beyond the horizon then all points of that string are assumed 
to be beyond the horizon. 

For the points that are visible an angle T is calculated 
according to the following formula: 


T = 


Tan -1 [ 


Sin (R) i 
S - Cos (R) J 


(2) 


0 

A polar angle about the subsatellite points is calculated 
as follows: 


U = 


Tan 1 [ 


Cos (YL) Sin (XL) 


Cos(LTSS) Sin(YL) - Sin(LTSS) 


] 


COS (YL) COS ( XL) 
(3) 


The Calcomp plotter uses Cartesian coordinates so these 
polar co-ordinates must be changed to the plotter coordinates, 
XP, YP, as follows: 


XP = T* Sin (U) 

YP = T*Cos(U) (4) 

An additional feature of this program is that it will 
draw the contours of constant elevation on the plot for any 
set of input elevations. 

The only inputs necessary for this program other than 
the data tape are the subsatellite longitude and latitude, 
the magnifying factor and the set of elevations for the con- 
stant elevation contours. Figure 6 shows a plot of this type 
with 70°, 45° and 5° elevation contours. Figure 7 is a 
flowchart of the program. 

The program can handle any number of cases in one run. 
The only limitation is the amount of CPU time specified on 
the job card. On the IBM 360/50 system a single case takes 
approximately 45 to 55 seconds, depending on subsatellite 
location and number of elevation contours to be drawn. 

A listing of the program is included in Appendix D. 
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Flowchart for the Polar PERSPECTive Program. 


Figure 7 
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VI. FOOTPRINTS OF ATS-F S-BAND ETV OFF-AXIS DUAL FEEDS 

NASA's Applications Technology Satellite-F (ATS-F), 
scheduled for launch in April, 1974, is going to test the 
feasibility of the direct delivery of ETV signals to low 
cost terminals in the newly allocated 2500-2690 MHz band 
in Alaska, Appalachian region, and the Rocky Mountain States. 
The antenna coverage program, described earlier in Section 
IV, has been modified to plot the S-Band footprints produced 
by the ATS-F spacecraft. 

The basic ATS-F spacecraft includes a prime focus feed 

( 2 ) 

complex having crossed-feed elements. The two ETV beams 

are generated by the 30 foot dish onboard the spacecraft from 
feeds which tie on the satellite North-South axis. Neither 
feed lies on the antenna boresight, the separation between 
the boresight and the nearest feed being about 0.7 degree 
and the beam separation being 1 degree. Since the beamwidth 
is also approximately one degree, the total coverage at any 
one time consists of two footprints lying in approximately 
a North-South relationship, the exact arrangement depending 
upon the subsatellite point and the boresight location. 

These factors have been programmed into the antenna 
coverage program described in Section IV. The inputs to 
the ATS-F S-Band ETV dual-feed coverage program (ATSFS) are 
the subsatellite and boresight locations, beamwidth dimen- 
sions and the plotting parameters. The plotting parameters 
are FX, FY, DX, DY and FACT, all of which are described above. 

The modified coverage program can be used to plot foot- 
prints either over the WORLDMAP or the MINMAP. Figure 8 
shows the footprints with a subsatellite location of 94°West 
and boresight pointed towards 81°West and 37°North. A 
listing of the program is included in Appendix E. 
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Figure 8. ATS-F S-Band ETV Footprints. 

Satellite at 94°W Longitude) 
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Appendix A 
WORLDMAP Program -1 

C WORLDS AP 

C 

0 I MENS ICN CCGRD UC24 ) » I BUF ( 1COO) ,PLTX ( 5 14 ) » FLTY (5 14 ), N< 23) 
Call PLOTS! IBLF, 1DOO ) 

C Rt AC SCALE FACTOR 

REAC(5,1CC) FACT 
iOJ FORMAT { F 1 0 . 3 ) 

CALL FACTCMFACTJ 
CALL PLCT( C .0,-29. 5 ,-3 ) 

F VX -- 3 . 1415926 

FVY=FVX/2.0 

C A X = 2 .*3 .1415926/50. 

C AY=3 . 1415926/26 . 

FX=-18C. 

F > =-9 C . 

GX=36C./50. 

CY=180./26. 

CALL PLCT ( 2 .0,2 .5, 23) 

C CixAki AXES 

C/LL AX I${ C .,0 . , 'LONGITUDE • ,-9 ,50 . , C . , F X , DX ) 

CALL AX1S<C.C,0.0, 'LATITUDE. * , 8 , 26 . , 90 . , F Y , D V ) 

CALL A X I S( C • 0 » 2 a • , 'LONGITUDE * ,9 ,50. ,C. , F X , D X ) 

CALL A X I S ( 5 0 . , C . , ' LATITUDE ' 0 ,26. , 90 . , F Y , 0 Y ) 

H c A C ( >3 ) N 
K = L 

Si 3 LAO COORDINATE BLOCK 
?J PE A C ( 8 ) CCORO 
8 i L = 1 
15 1=1 

IF ( COORD (2*L-1 ) . Gt .4 .0 )L = L + 1 
3 7 PLTX ( I ) = CGCP0(2 J ''L) 

PLT Y ( I ) = CCC Rl) ( 2 *L- 1 ) 

C CiifctK FOR END OF STRING 

I F C PLTX ( I ) . EQ.O. O.ANC.PLTY (!) . EQ.C .0) GG TC 47 
I F ( I • L T . 2 ) GC TC 31 

C C.luCK FOR -18C TC +180 kRAP ARfLNC 
TESl = A3S(PLTX(I)“Pl. TX(I-l) ) 
l F ( TE ST .GT .2. )G0 TC 27 
31 L = L + 1 

1 = 1 + 1 

OC 10 37 
Hi PLTX ( I ) = F V X 
PLT V ( I ) =FVY 
PLTXU + 1)=CAX 
PLTY ( 1+ 1) =C/.Y 
NPT S = I-1 

C OKA* STRING OF PCINTS 

CALL L ! NE ( PLT X, PITY, NPT S» 1.0, C ) 

50 L=L+1 

. 1 F { I. . L E • N ( K ) ) GG TG 15 

X = K H 

IFtK.NE. 12) GO TO 5 
CALL PLCT (G. ,C . ,-3 ) 
j IF( K.GT.23JGC TO 90 
GC TU 7 C 
cl PLTX ( I ) = F V X 
PLTY ( L ) =FVY 
PLTX ( I + 1 ) =CAX 
PLT Y ( 1 + 1 ) =0A Y 



Appendix A 

WORLDMAP Program -2 


NPT S = I- l 

C U*Afc STRING CF PCINTS 

CALL LlNt(PLTXfPLTY»NPTS»ltflfC) 
G G TO 15 

5J CALL PLOT(0.«0.«-3) 

C CkAI* GRiC LINES 
UC 45 J= 1 » 2 5 
X = F LG A T (2*J-1 ) 

N RG * = J- ( J/ 2 ) #2 
IF(NRE V .£C.C) GU TO 46 
CALL PLCT(X,0.,2) 

CALL PLCT (X ,26. , 2) 

GC IG 45 

4o CALL PLOT (X, 26. ,2) 

CALL PLCT(X,0.,2) 

4d CONTINUE 

DC 45 J = l , H 
Y=FLC AT (2* J-I ) 

NREf = J- ( J/2 ) *2 
iF(NRG^.EC.O) GO TO 48 
CALL PLCT(C.,Y,3) 

CALL PLCT < 5 C . , V ,2) 

GC TO 45 

48 CALL PLCTI50. ,Y,3) 

CALL PLCT(0.,Y,2) 

4 ) CCNTINLfc 

CALL PLCT(C.,C.,595I 

STOP 

cNO 



Appendix B: 

MINMAP Program -1- 


C i v I N f“ A P 

C 

DIVENS [GN IUUF(IOCG) , CCD RD ( 1 02 4 ) , N ( 2 2 ) , X F ( 5 14 ) , YP 1 5 14 ) 
Call PtCIStIBUF.IOOC) 

CALL FLCTIO. ,-29.5,-? ) 

C Kb AC SCALE FAC T C R 

Rt«C {*5,9 9) FACT 
Si F C R N A T I F 1 0 . 3 ) 

CALL FACTCPCFACT) 

CALL PLCT(l.«.5»23) 

CCN 'vTR-3. 1415 S 265 3 5/100. 

C RE aD rJCLiN C AR I F S CF MAP 

REA0(5fiCC) XLL »XUL » V L L * Y U L 
iCJ FCRNAT 14F10.3) 

F VX = XL L 
F V Y = Y L L 

C A = ( VLL-VLL ) / 1 0 . 
rtRITE(6,20C) FVX,FVV,CA,FACT 

2 CD FCKf''AT{lX, , FVX=',FlC.4,4X,*FYY'=»,F10.4,4X,'CA=',F10.4, 
* 4 X » ' F A C T — * fFlO.A) 

X A X L = ( XUL-XLL ) /C f> 

V A X L = 1 0 . 

C OK Art AXES 

CALL AXIS (C. ,C. ♦ ' LCNG I TUOE ' ,-<3 ,XA XL , C . , F V X , C A ) 

CALL AXlS(C.»C.t 'LATITUDE , t8fYA.XL»9C«fFVV»DA) 

CALL AXIS(C.,Y/-XL, ' LCNd TLDfi ' , 9 , X A XL , C . , F VX ,C A ) 

Call A>lS(XAXL,0.f 'LATITUDE • ,-8,YAXL,SC.,FVY,CA.) 
XLL-XLL*CC NY T R 
C A=CA*CCNV TR 
XLL =XLL*CCNYTP 

V L L =YLL*CCNVTP 
YLL=YUL*CCNV1R 
RE«0(8) N 

DC 1C K = 1 » 2 3 
READ (8) CCCRD 
NL = 1 
15 N P= l 

IF(CQCRC(2*NL-l).LT.4.»rc TT 11 
N L = N L +■ 1 

li XF(NP)=CCCKD(2*NL> 

YF{NP)=CCCRC(2* M— 1 ) 

C ChbCK FCR END CF STPING 

I F C >P ( \P ) . EC .0 . . AND . YP ( NF ) . EC . 0 . i GC TC 12 


li-UCK FUR PCiNT IN BCUNCS 


IF(XP(NP).LT.XLL)GC 

TO 

1? 

IF(XP(NP).CT.XLL)GC 

TO 

1 3 

iF(YPINP). LT.YLL )GC 

TO 

12 

IF( YPCNP) .CT .YLL )GC 

TG 

1 3 

N L = N L + 1 



N F = N P + 1 



IFIM.LE.MK) ICC TC 

l 1 



GC TC 1C 
12 XF(NP)=XLL 
YP ( NP ) - YLL 
XP { NP+1 ) =D A 
YP ( NP+ 1 ) =L A 
N FI S=NP— 1 

C DKArt STRING CF FCINTS 

IF (NP7S.LT .2)GC TG 14 



Appendix B 
MINMAP Program -2 


CALL LINEUP, YP,NPTS,1,C,0) 
l* f\L = NL + l 

IF( NL .GT.N 1 K) >GO Tf 10 
GC TO 15 
13 XF ( N P ) - X LL 
YP(NP)-YLL 
XFINP+l)=DA 
YF ( NP+1 ) = DA 
NF7 S = N P - 1 

C STRING CF FCIN7S 

lFINPTS.LT. 2 I GC TO 17 
CALL L I NE ( XP,YP, NPTS , 1 ,0 ,0) 

17 N L = N L + 1 

X = CCCRC ( 2*N L ) 

Y=CCCRD(2*NL-1 ) 

C C.lt-U K FCR LND CF STRING 

IF( X.E5-.C. .ANC.Y.EC.C.IGC TO 1 A 
C CHECK FCR PCINT IN PCUNCS 
1FU.LT.XLL1GC TC 16 
IFO.GT.XUDGC TG 16 
IFl Y.LT.YLDGC TC 16 
if ( Y.GT.YUL1GC TC 16 
NP= 1 

GC TO 11 

lo IF1NL.LT. MK) ICO TC 17 
10 CCNTINUE 
C OKAW GRIC LINES 
DC 2 0 1 = 1, 9 
Y =F L C AT ( I ) 

CALL FLCT1C . ,Y ,3 J 
20 CALL FLCT { X A X L , Y ,2) 

N > = I N T 1 X A X L ) 

IF1FL0AT (NX) .EC.XAXL ) N X = N X— l 
DC 30 J= 1 , N X 
X=FLCAT( J) 

CALI PLC7U,0.,3> 

3 J CALL PLCT(Xtl0.t2) 

CALL FLCT(0.,0.,999 > 

STOP 

END 



c 

c 


20:? 


2Gi. 


2 

9CJ 

i 

j.C J 
5Qi 


ANTCVL V 


Appendix C 

ANTOVLY Program -1- 


I^PLiCIT REAL*S(A-Z) 

R <i A L * 8 C S i N » DCC S,DS£RT» 
I NT EGER *4 NSTEPS » ICT » I» 
RE A L *4 X P ( 5 1 4 ) * Y P ( 5 1 4 ) 

1 NTcG e8*4 J ,Jl*L»K f NFL 
INTEGER *4 FM(4 )/'(', ' 
INTEGLRX-4 C I G I T ( 10 ) / • 1 ' 


CARS IN ,DARCCS,CATAN 

N.NPLTf IRUM4CC0 ) ,NP 1, NP2 . NCP ,LOR 


• , 'F7.3* , ' ) * / 
» , 2 , t , 3', , 4 , » , 5'» , 6 , ♦ ' 7 ' 



RL3L*4 LCN( 722 ),LAT{72?),M,XCT,YCT,FX,FY f DX,DY,XSS,YSS 


' 10 •/ 


INTEGER I M , I P * N P 
INTEGER#4 1 S YN ( 5 ) 

CATS ISYy/ 3 , 4 , 1 L , 9 , 10 / 

DINENSICN CE ( 10) ,DRS( 8) ,RRWS(8 ) 

DATA C3S/.1».2,.5,1.*1.5,3.,5.,1C./ 

LATA RPWS/. 18 , . 26 , . 4 , . 56 ,. 7 , 1 . , 1 . 27 , 1 . 7 / 

DATA PI/3. 141592(53 5E5 793DO/,RE/3 .96C3/,C IS1/2.6 2604/ 


N = C 

RtA 0 ( 5,205) FX»FY,cX,2Y,rX,CY 
FCRPAT (6F1C.3 ) 

CALL PLCTS ( IBUF ,4000 ) 

8 1 A C l 5 * 20 1 ) FACT 
FCR«AT(F10.3) 

CALL FACTCRt FACT) 

CCN\iTR = PI/180. 
i3M=C APS INC RE/C 1ST) 

CHK =81 . 3*CC NVTR 

HuAC ( 5,9CC ,£ND= 1000 ) NPLT 

FCRNAT (II ) 

DC 23 NPL= 1 , NPL T 

READ! 5 ,5CC , £ ND = 5 5 0 ) CL C N S S t DLCNCT , D L AT C T ,DINCR,OCELT 
FCRNAT (5F9.3 ) 

RE4C(5,5Cl)TFETAl,BWlA,RV»2A f L 
FCRNmT (3F9. 3 , 12 ) 

F IY I ( 2 ) - C I G IKLI 

RE A0( 5 , FNT )( CB ( l ) , I = 1 ,L ) 


NC 8- C 


N = N*1 

WRIlE(6,6CC)N,CLCNSS,niCNCT,CLATCT,CCELT,PW2A,BKlA,THLTAl 
oOj FORMAT {// • 10ATA SET NO. ',I2/4X,*SUR SAT LONG ='*F9.3/ 

1 -,X,'3CR SGHT LCNG = ' , F 9 . 3 / 4 X , • B C R SGFTLAT =',F9.3/ 

2 4X» •OtCLlN/.TICN = , ,F9.3/4X,'MN BV fcC TH • , 5 X , • = • , F 9 . 3/ 
i 4X,»N*X BVWDTHS5X, •= SF9.3/4X,' ORIENT ATI CN =',F9.3) 

Q CCiWERI FRON DEGREES TC RADIANS 
INCR=OINCR»CCNV TR 
LCNCTRMCLCNCT-DLCNSS )*CCNVTP 
LATCTR=CLATCT*CCNVTR 
DELI=CCELT*CCNVTR 
Bfcl=BWlA*CCNVTP/2.0 
Bk2 = DW2A*CCN\/TR /2 .0 
T F E I A = T H E T A i*CCNVTR 

C CHECK FCR BORE SIGHT LCCATICN IN RANGE CF SATELLITE 
AL = CAPCCS (DCCS (LONG TP >*CCCS ( I. A TC TR + DFLT ) ) 

IFl AL.LE.CFK )GC TC 21 
20 WRITt(6,4GC) 

•*GJ FCRIYF. T 1 1 H 0 * 'BGRESIGFT LCCATICN IS NCT IN RANGE CF S AT FLL I TF . * / ) 
GC TO l 

C CUMPLTE VECTCR FRCN SATELLITE TC PORES IGFT LOCATION 
kSX=RE*CCCS(LATCTR)#CSIN(LCNCTR) 



Appendix- C 

ANTOVLY Program -2- 

RSY=OIST*CCOS(CELT)-PE*CCCS ( L ATCTR )*CCCS (LONC TR ) 
*SZ=OIST*OSIN( CELT >+RE*OSIN(L ATCTR} 
RSP-USGRT(PSX**2+RSY**2+RSZ4<2) 

GAM=QAPSIN(PE*DSIN(AL)/PSM) 

EL=((PI/2.0)— (AL*GAR))/CCNVTP 

WRITE(6,300) EL 

jCJ FCRI''AT(AX» , ELEVATICN , »6Xf , = , *F9.3) 

G£NCM=OSCRT (RSX<'*2+PSY 5 ! [ ^2 ) 

P1TCH=CAPSIN (RSX/DENC^) 

HCLL=OATAN ( fiSZ/DENCN ) 

G P = OC GS (PITCHI 
CP = CCCS (PCLL ) 

SP = QSIMPITCH> 

S R=C $ I M RC L L ) 
ii DC 32 K = 1 » L 
NCB=NCS+1 
LC3=L-K* 1 

C INTERRELATE FCP RELATIVE RFANVuICTF 
DC 8 J=1 » 8 

IF{CE(LCe)-CES(J)>6t7,3 
/ FEVi=RBWS(J) 

GC IG c - 
jCCMiNLc 

RBw = (ce(Lcei-DBS(Ji))>MRBws(j)-PBus(Ji))/(res(j)-ces(Ji)) 

i +R3RSIJ1) 

■i A-RSM*GTAN(BV«1<RPW) 

B=RSM*0TAN ( EW2*RBW ) 

ALPHi = !3WlA*RRW 
ALFH2= PR2A*«ER 
NSTEPS=(360.3/CINCR)+1 
BRWC TH = CALPHA’>P3Vs 
OC IC l C T = I ,NSTEPS 
3 E T A N = ( ICT-l)*INCR 
A N G = tJ c T A N + T F E T A 

UM= l.O/CSGRT ( ( CCCS l BET AN )/A )**2 + ICS !MRETAN)/B)**<) 

C CJMPLTE N-TH R VECTCR FRCH SATELLITE TC LCCLS CN EARTF 

MNPX = RSM*CR*SP + UR*CCCS ( ANG >*CP-UR*DSIN ( ANG) *SP*SP 
— MfvPY= f JS^«CP ,> CP + UN*CCCS(ANG) <! SP-UY*CSIN(ANG)*SR*CP 

MNPZ = RSN’ )t SB + LI' / *OSIN(ANC:)*CR 
RNX-MNPX 

FNY=HNPY*CCCS(CELT )4PNPZ*CSIN (CELT ) 
f“'NZ = MNFY ! M-CSIN(CELT)) + RNPZ*CCCS(DELT) 

8N-CARCOS (MNY/CSCRT (VNX**2 + RNY* ,! ‘2 + PNZ* ,! '2 ) ) 

I F ( EN.GT .BR > GC TC 3C 
CN=Pl-CARSIMDSIN(BM*D!ST/RE) 

DN=P l— (BN+CN ) 

HNL=OSGRT(RE**2+LIST**2-2.*RF*OTST*CCrS (CN) ) 

GENCi v l = DSCRT (MNFX**2+RNPY**2 ) 

P 1 1 ChN=OAk 3 IM f^NPX/CFNCM ) 

RCLLN = P. ATAN (MNFZ/DENCM ) 

c CJrBPLIE N-TH VECTOR FPCN CENTER OF EARTH TC LOCUS 
REl=MNL*CCCS(RCLLN)*OSIN(PITCHN) 

KEJ=DI ST*CCCS( CELT l-MM+CCCS (PC LI N ) *DCC S(PITCFN) 

Rfc:K=RNL*CS IN ( RCLLN )-PI ST*OS IM Ot LT ) 

GC TG AO 

q If M-TH V V EC 1 C R DCFS NOT INTERSECT FARTF CCYPLTE VECTOR 

c FkCr center cf earth tc forizcn seen ry satellitf 

3D T AUN-CARCCS (RNZ/CSQRT { l«NX** 2 + R N7-** 2 ) ) 



Appendix C 

ANTOVLY Program -3- 


IF < I* NX .LT.O.O)TAUN=-TAUN 
KEl-RE*nCGS lev ) *DS I N ( TAUN) 

FEJ=RE*(-CC.CS(CLLT )*CSJM6M )-T S I N ( CELT ) * ECO S< P* l*CCCS( T*UN) > 
REK=kE*{-CS IN( CELT )*DSIN(BM I+CCCS {0ELT)*GCOS(9R )*CCCS (TAUN) ) 
ttj CEN=DSCPT(PEI**2+REJ**2> 
r CuMPLTE L C NG I TLCE mNC LATITUDE CCCRCIMTtS 
LAT (I CT ) -SNCL (DATANIRFK/DEN I/CCNVTP > 

10 LCM ICT ) =SNGU CAPS INIRF I / CEM/CCNVTR + CLCNSS ) 
wRI TC (6 .ICO) L3 ( LCB ) ,M.PF1 ,ALPH2 

iGJ FCRRATI1H— »3X.'AT ‘.FA.l,* TO L E V E L : • / A X , • M AX RPWCTM '»F5.2/ 

L AXt'RIN BRkOTF= • » F 5 . 2/ / 

Z 5X,‘LCNGITUDE',7X,*LATTTUrL'/5X.' (CECREESJ'.eX,' (CEGRFES)' 

3 /) 

DC 80 I=1,NSTEFS 
URI TE ( 6 » 2 C C ) LCN( I) tLAT ( I> 
zCJ FC^aTUH ,4X,F9.3.6X,F9.3) 

IF (LCN( I > .LT .-18C. ) L C M I )=360.+t.GN( I ) 

IF (LCN ( I ) . LT .F X ) LC M I) = FX 
iF(LCNtn.GT.ex) LCMT) = FX 
i F ( l AT ( n.LT.FV) LAKD-FY 
IF lLAT { I ) .GT .E Y ) LAT{ I ) = EY 
SuCC'^TiNUE 

CCiWT R = S N G L ICG fv V T R ) 

XCT=( SNGL( CLCNCT J-FX )/HX 
YCT = ( SNGLICLATCT )— FY)/CY 

CALL SYRHCL(XCT,YCT,0.14,ISYR(N>,C.,-1) 

C i*4_ C T CCCRCIFAIES CF LCCLS 

1 F = 1 
5j IP=1 

11 I F ( ABS (LCM 11* + l )-L CIV ( I V ) ) .GT . IOC . )GC TC 70 
3/ XP (1 P )=LCM 1M» 

YP( IP ) =L AT ( IP ) 

IMIM + l 
IP=IP+1 

I F ( IF. LE. NS 1 EPS) GO TC 47 
X P ( I P ) = F X 
XP{ IP + 1 ) -D X 
YP( IP) 

Y P ( IP+1 ) =D V 
NP* IP- 1 

CALL LiryEIXP.YF.NP. 1.0,0) 

GC TG 32 

4/ i F ( IM.ti. N STEPS ) C-G TC 37 
GC 1C 31 

70 XP( IP ) = LCM IM 
YPIIP) = LATUM) 

XP ( IP + 1 ) = F X 
XF{ IP + 2 ) = DX 

Y P ( IP+1 ) = F Y 

Y P ( IP + 2 ) = C Y 

CALL LINEIXP.YP, IP. 1,0,0) 

IR=IR+1 

GC TQ 50 
8Z CCMIME 

XSS=(SISGL(CLCNSS)-FX)/CX 
YSS = ( SNGL (-CDELT )— FY )/CY 
2.0 CALL SY, v BCL{XES.YSS»0.21,ISYV(N),C.,-l) 

0 5o CALL PLC T ( 55 . , — 1 • » 23 ) 

GC 10 2 



Appendix C 

ANTOVLY Program -4- 


C CLOSt PLCTTAFE 
LcCJ CALL PLCT (C. ,C. ,<599 ) 
rfRI T ^ (6,601 > 

oCi FCPPAT ( 'CPLCTTAPE CLCSEO*) 
S7CP 
EN u) 



Appendix D 
PERSPECT Program 


C FE R SP L CT 

C 

DIMENSION I OUF { 1COO 1 ,CCORD ( 10 24 ) ,N ( 2 ; ) , X F ( 5 14 ) , VP C 5 14 ) 
.DIPcNSICN £ l ( 1 0 ) 

REAL*'-, LNS S » L T S S 
PEAL*4 LNSSR 

CM S PI ,S/3. 1415026,6.6166/ 

CALL FLCTS ( IBUF , 100C ) 

CALL FLCTIC.,-29.5,-2) 

R. c A C ( 5 >400) FACT 
4C0 FCRMTIF10.3) 

CALL FACTOR! FACT) 

CALL FLC7 ( 17. , 15. ,23 ) 

CC'NVt TB = P I / 180. 

TIM>=ARSIM 1 . / S ) 

F V = C • C 
C A = CC N V T R 

ANAX=SC. C N V T R — TMA X 
C RtAC SUeSATtLLlTL CCCRHI NATES 

j READ! 5 ,1CO,END=ICOO)LNSS ,1 TSS 
100 FORMAT (2FLC.2) 

LNSSP=LNSS’* t CONVTR 
CELTR=LTSS *CC N V T P 
SINO-=S INIOELTR ) 

CCSC=CCS(UtLTR > 

REAC131N 
K = 1 

£0 P CAC ( S ) COORD 
L = 1 

10 1 = 1 

L F ( CGCRO (2*L-i ) .GE.4. )L=L+1 
YL=CCCR0(2*L— 1 ) 

XL = CCC»D(2*U-LNSSP 
SINLT = SIN(YL) 

CCSIT=CCS( VL ) 

SINLN=SIN( XL) 

CCSLN=CCS( XU 

CC$R=$IND*$INLT+CCSC*CCSIT*CCSLN 

R=ARCCS (COSP ) 

DIF=RRAX-R 

C IS FIRST PCINT IN STRING V I TP I N 30 CECREFS CF HCPIZCN? 

I F ( ABS(CIF) ,LT. .52359881GC TC 20 
C FIRST PCINT >30 CfcGREES F RC R FCRIZCN... 

..C . ... IS IT VISIBLE? 

IFICIF.GT.C.IGC TC 70 
4u L = L+ 1 

IF ( CCQRC (2*L-1 ) .FC.C . .ANC.CCCRC ( 2*L ) .FC .0. ) CO TC 60 
GC TC 40 

C FIRST PCINT <3C CECREES FRCR FCRIZCN... 

C_ ... CHECK EACH FCINT FCP VISIBILITY. 

20 I F ( C I F . G T . C . )GG TO 70 
C NCT VISIBLE.. .CHECK NEXT PCINT 

L = L ♦ 1 

IFICC0RDI2AL-1 ) .EG.0..AN0.CCCRC(2*L).EC.C.)C0 TC 60 
YL=CCCRD{2*L-1 ) 

X L = CCC RC (2*L)-LNSSR 
5 I N L T = S IN ( YL ) 

CCSLT =CCS( YU 
SINLN=S I N ( X L ) 



Appendix D 
PERSPECT Program 


CC S LN = CC S ( X L ) 

CCSR=S IND*S INLT*CCSC*CCSLT*CCSLN 
R = A FCC S I CO S F ) 

C1F = HMX-R 
GC 1C 20 

C VISIBLE. ..ChtCK NeXT FCINT 

7 J T=ATAMSIN(R)/(S-CC$R)) 

L=ATAN?(COSLT*SI N L N , C C SO* S I M T-S I NC*CC S l T *C CS L N ) 
>PlI) = T*SiMU) 

v f ( n = r*ccs (u> 
i = i*i 

L = L* L 

IF (CCC«D(2*L-1 ) . EQ .0. .ANC.CCCRC (2*L ) .EG .C.) CO TC 50 

YI = CCCRD(2*l-1 ) 

>L=CCCPD (2*L)-LNSSR 
SIMT=SIMYL) 

CCStT=CCS(YL) 

SINLN=SIN(XL) 

CCSLN=CCS( XL ) 

CC SR=SIND*SINLT+CCSC*CCSLT*CCSLN 
k=ARCCS (COSP) 

CIF = RMX-R 

If (C IF.GT.O.JGC TC 7C 
L -L - 1 

5J XP ( I) =FV 
> P ( 1+ 1 ) = DA 

Y F ( l ) = F V 

Y P ( I ♦ 1 ) =CA 
N P T S = I - 1 

C PLC1 STRING 

CALL LINE(XP»YP,NPTS,1,C,C) 

6 J L = L ♦ 1 

lF(L.Lc.N(K ) )GC TO 1C 
K =K + l 

I F ( K .LE.23 )GG TO 9C 
C GRAlft HCRIZCN 

OC 50 J = 1 , 7 3 
GET A = 5 • * ( J— 1 ) * CC NV T F 
XP( J)=TPAX*CCS(8FTA) 

5o YP( 2 ) =TNAX*S IN ( BETA ) 

XP( J+l )=FV 
X P ( J-* 2 ) = DA 
YP(»j + 1} = FV 

Y P ( J + 2 ) =DA 

CALL LINEUP, YP, 73,1, 0,0) 

CALL STPBCLtO. ,0., .2B,3,C.,-1 ) 

C OPAw CONSTANT E LEV A T I C N CCNTCIJRS 

FlAQ( 5,200) NTS 
200 FCRNAT(Il) 

IF ( NT S.EC. 0) GC TC FA 
REA3(5»300) (cL( I )» I=1»NTS) 
j Q j F C R w AT (9F5 . 3 ) 

OC 5 t K = 1 , N T S 

ANG2= 190. + ELIK) )*CCNVTR 

AN G 1 = A R S I N l SIM ANG? )/S) 

OC 55 2=1,73 
6ETA=5.*(J-1 ) * C C N V T F 
X P ( J ) = ANG1 *CCS { BET A ) 

9 j YF(j)=ANGl*SlN(BfTA) 



Appendix D 

PERSPECT Program -3- 


X P ( J+ 1 ) =FV 
XP( J + 2 ) = Da 
YP(J+1)=FV 
Y P ( J + 2 ) =C A 

CALL LINE(XP,YF,73,1,C,0) 

YN=ANGl/DA+.05 

CALL NUPBER(0.,YN,.21,EL(K),0.,-1) 

CCMINIE 

CALL SYPBCL (-1 .33,-10. , .2 1 , 'LONGITUDE = ' . 0 . . 1 1 ) 
CALL M M 8 1 R (ASS. ,999. ».21»LN$S»C. »?> 

CALL SYMBGL(— I. 33, -10. 35..21,' LATITUDE ='»0.,11) 
CALL MMBtR (909. .90S. ..2 1.LT5S.0. ,3) 

CALL FLCT ( 30. ,C . ,-3 ) 

REWIND 8 
GC IU 5 

it C J CALL PLCT (20. ,-15. ,990 ) 

STOP 
c NO 



Appendix E 

ATS-F Dual-Feed Program -1- 


C 

c 


ATS-FS 


99 


IMPLICIT REALMS! A-Z) 

REALMS DSIN»DCGScDSGRT,CAR$IN*OARCOS»CATAN 

I NT EGER $4 NSTEPSi?ICT*I*N»NPLT*IBUFtlOOO)»NPi»NP2*NC8*LD8 

RE AL*4 FX*FY*DX*OY 


INTEGERS 

INTEGERS FNTf^/M*** * * * F 7 . 3 * » * ) * / 

INTEGERS 0IGI¥5i05/*l , »'2'» , 3 , **4 , » , 5 , »*6* »' 7* » *6* 
REAL $4 LON (722 I * LATf 722 ) »N 1 * XCT *YCT 
INTEGER U 


' 9 ' * 


REALMS ROLL! « ANGDIF 12 I 
DIMENSION OBSIO) *DBS(8)*RBMS(8I 
CATA C8S/olj *29 »5 * l.*1.5*3. *5.*1C./ 

DATA R8WS/® 18 » . 26 * .4 * . 56 * . 7*1 . * 1 . 27 * 1 . 7 / 

DATA PI/ 3. 1415926535 897 93 DO /*RE/3. 96 C3/»0IST/2.62 604/ 

R E A 0 ( 5 *99 } FACT 
FORMAT ( F I0®3 J 


10 */ 


N = 0 

NPLT*0 

CC N VT R = P I / 1 80 • 

ANGOIFCl )=.5^CCNVTR 
ANGC KFI2)*-.5*CCNVTP 
BM=DARSIN(RE/OISTI 
CHK-81®3*CCNVTR 

1 READ ( 5*500 *ENO= 1000 )DLCNSS ♦ DLCNCT * CLATCT »D I NCR » DOELT 

500 FC RN A1 ( 5F9 ®3 1 

PEA D I 5 * 501 3 0W1A,BW2A 

501 FORM A? 5 2 F9 ® 3 1 

WRITE (6*600) N,OLCNSS,OLCNCT,CLATCT * DDEL T . BW2A * BW 1 A , TH ET A 1 
600 FORMAT C //• 1DATA SET NO. »*12/4X*'SUB SAT LONG =*,F9.3/ 

1 4 X* 0 BOR SGHT LCNG = ' ,F9. 3/4X* ' BOR SGHT LAT =*,F9.3/ 

2 4X» "DSC LI NATION = • , F9. 3/4 X, • M IN RMViCTH • • 5X , ' * • ,F9 . 3/ 

3 4X, 8 MAX BMhDTH',5X**=*,F9.3/4X, 'ORIENTATION *',F9.3) 

C CCNVEPT FROM DECREES TC RADIANS 

INCRaCINCR^CCNWTR 

LCNCTR=«DLCNCT-OLCNSS)*CCNVTR 

LATCTR=CLATCT*CCNVTR 

OELT^CCELT #CCN VTR 

BU’BWIA+CCNVTR/a.O 

eV2=BW2A*CCNVTR/2.0 

C CHECK FOR BORES IGHT LOCATION IN RANGE OF SATELLITE 
AL=OARCQS(CCCS<LCNCTP)*CCOS (LAT CTR+DELT ) ) 

IF( AL.LE.CHKl GO TO 21 

20 W R I T E ! fc * 40 C 5 „ . .. 

400 FORMAT ( 1H0 * 8 BOR ESI GHT LOCATION IS NOT IN RANGE OF SATELLITE.*/) 

GO TC 1 

C CQMPITE VECTOR FRCM SATELLITE TO BORESIGFT LOCATION 

21 RSX«RE*0CQS(LA1CTR )*DS IN(LONCTR) 
RSY=DIST*DCCS(DELT)-PE*CCO$(LATCTR)*CCCS(L0NCTR) 
RSZ=OIST*DS!N(OELT S*RE*DSIN (LATCTR) 
RSM=0$QRT(RSX**2*RSY**2+RSZ**2> 

GAM=CAR$IN!RE$DSINCAL)/PSN) 

El-f ?PI/2.0)-( AL-&GAM) l/CCNVTR 
MRI1 E(6 «300B EL 

300 FCRWA7 84X. 'ELEVATION* *6X»*- »*F 9. 3) 

DeNCM=0S0RTSRSX**2*>RSV**2 ) 

PITCF=CARS1N( RSX/OENCM) 



Appendix E 

ATS-F Dual-Feed Program -2 


RCL L 1 = DATAN ( RSZ/CENCM ) 

OC 31 IJ=1,2 
RCLL=RCLL1*ANGDIFC I J ) 

CP=CCC S ( P I TCH ) 

CR = CCGS< ROLL) 

SP 3 CS IN I PITCH! 

SR = CS I M ROLL ) 

NPU = NPLT+1 
IFCNPLT.GT.11GO TO 15 

C FIRST PICT, OPEN PLCTT APE AND SET ORIGIN 
CALL FACTOR (FACT ) 

CALL PLOTSUBUF, 1000) 

15 CONTINUE 

A=R!M*OTANI BW1 ) 

B=RSM*CTAN(8N2) 

A L PH 1=8 W 1 A 
ALPH2=BW2A 
UM= A 

I F { I J . EQ . 2 ) UM = B 
N STEP S= 5 360.0/DI NCR )*1 
DC 10 ICT-1 V NSTEPS 
BETAN=( ICT— 11*1 NCR 
ANG-BETAN 

C COMPETE N-TH M VECTOR FROM SATELLITE TO LCCliS CN EARTH 
MNPX=RSM*CR*SP+UM*OCOS( ANG)*CP-UM*DSIN( ANG)*SR*SP 
MNPY=RSM*CR*CP+UM*DCOS< ANG1*SP-UMCSIN ( ANG)*SR*CP 
MNPZ=RSM*SR+UM*DSIN(ANG)*CR 
MNX=MMPX 

MNY=MNPY*OCCS(CELT)*MNPZ*OSlN(CELT) 

MNZ=MNPY*I-OS INI DELT) ) +MNPZ*CCOS (CELT ) 
8N=CARCOSiMNY/OSQRT(MNXA*2+MNY**2*MNZ**2) ) 

IF(BN.GT.BM)GO TO 30 
CN=FI-DARSIN(CSIN(BN)*DIST/RE) 

DN= P I- ( BN+CN) 

MNL=0SCRT<RE**2+DIST**2-2.*RE*DIST*CCCS(CN) I 
OENCM=OS<3RT{MNPX**2 + MNPY**2 ) 

PITCHN=OARSIN(MNPX/DENCM> 

RCLLN=DATAN(MNPZ/DENCM) 

C COMPETE N-TH VECTOR FRCM CENTER OF EARTH TC LOCUS 
R£I=MNL*CCCS4 RCLLN l*OS IN CPI TCHN I 
REJ=D I ST*DCOS ( DELT )-MN L*DCOS (PC LLN )*CCOS 4P I TCHN ) 
REK=MNL*OSIN4RCLLN)-DIST*DSIN(OELT) 

GC TC 40 

C IF N-TH M VECTOR DOES NOT INTERSECT EARTH COMPUTE VECTOR 
C FRC M CENTER OF EARTH TO HORIZON SEEN BY SATELLITE 
2d TALN=CARCOS4MNZ/CSQRT(MNX**2+MNZ**2l ) 
IF4MNX.LT.0.0)TAUN=-TAUN 
REI=RE*DC0S4BM)*DSIN(TAUN) 

REJ=RE*(~OCOS(DELT )*0S I Nt BM )-DS IN I DELT ) *CCOS( BM )*CCCS4TAUN) ) 
REK=RE*4— OS INI OE LT )*OSIN(BM)*CCOS(DELT)*CCOS(BM )*CC0S4TAUN) ) 
40 0EN=DSCRT(REI*+2+REJ**2 ) 

C CCMPLTE LCNGITUOE AND LATITUDE COORDINATES 
LAT(ICT)=SNGL(0ATAN(PEK/DEN1/CCNVTR) 

10 lCNUCTi=SNGL(CARSIN{REI/OEN)/CCNVTR+CLCNSS ) 

CC eo I=l,NSTEFS 

m kRIl£I€,20CJ LC N( I I , LAT ( I ) 

200 FORMAT 1 1H , 4X , FS .3 , 6X , F<3 «3 ) 

FX=- 151 . 

FY = C. 



Appendix E 

ATS-F Dual-Feed Program -3 


0 = 4 . 

DY=4. 

NP1=NSTEPS*1 
NP2 = N STEPS+2 
LCMNPl)=FX 
LCN (NP2)=0X 
LATCNP1)=FY 
L AT ( NP2 l=DY 

C PLOT COORDINATES CF LOCUS 

CALL LINE(LCN*LAT*NSTEPS»1»0*0 ) 

21 CONTINUE 

CALI PL0K35. 23) 

GO TO 1 

10C0 IFCNPLT.LE.C) GO TO 1001 
C CLOSE PLOTTAPE 

CALL PLCTIO.O f G.C,999I 
1001 hRITE(6.700) NPLT 
7C0 FORPATt • 1NUHBER OF PLOTS PRCOUCEC =*»I-) 

STOP 

END 





_ INPUT FORMATS 


1 . 

ANTOVLY 



1 . 

FX, FY, FY, DY 

4F10.3 


2. 

Magnification Factor 

F10.3 


3. 

Number of beams 

11 



4. SS long, BS long, BS lat, 
Circumference increment, 
Declination 

5F9.3 



5. Orientation, BW1, BW2 , 
# of DB levels 

3F9.3, 12 



6. DB levels 

nF7 . 3 



repeat 4,5 & 6 for each beam 


2. 

MINMAP 



1 . 

Magnification Factor 

F10.3 


2. 

XLL, XUL , YLL, DUL 

4F10.3 

3. 

WORLDMAP 



1 . 

Magnification Factor 

F10.3 

4. 

PERSPECT 



1 . 

Magnification Factor 

F10.3 


2. 

SS long, SS lat 

F10.3 


3. 

# of constant elevation 
contours 

11 


4. 

Const, elev. contours 

nF5 . 3 


repeat 2,3 & 4 for each case 


5. 

ATSFS 



1 . 

Magnification Factor 

F10.3 


2. 

SS long, BS long, BS lat. 
Circumference increment, 
Declination 

5F9.3 


3. 

BW1, BW2 

2F9.3 



